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Abstract. We investigate how thermahzation of gluons depends on the initial 
conditions assumed in ultrarelativistic heavy ion collisions at RHIC. The study 
is based on simulations employing the pQCD inspired parton cascade solving 
the Boltzmann equation for gluons. We consider independently produced mini- 
jets with pt > Po = 1-3 ~ 2.0 GeV and a color glass condensate as possible 
initial conditions for the freed gluons. It turns out that full kinetic equilibrium 
is achieved slightly sooner in denser system and its timescale tends to satu- 
rate. Compared with the kinetic equilibration we find a stronger dependence 
of chemical equilibration on the initial conditions. 
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1. Introduction 

Recently we have studied kinetic and chemical equilibration of gluons in a central 
heavy ion collision at RHIC energy employing a microscopical transport model in- 
cluding multiple gg <-> ggg scatterings [^j. The results showed that even for an 
initially dilute system, dNg/dy ~ 200, when the initial conditions are assumed as in- 
dependent minijets with px > po ^ 2.0 GeV, overall kinetic equilibrium is achieved 
at 2 fm/c and full chemical equilibirum follows later at about 3 fm/c. In addition, 
the expansion of the considered parton system behaves (quasi)hydrodynamically. 
Especially the total transverse energy per unit rapidity at midrapidity, dEx/dyly^Q, 
decreases from initial 500 GeV to 270 GeV at final time 4 fm/c due to the longitu- 
dinal work done by the pressure built up during the thermal equilibration. 

Since the experiments at RHIC at y/s = 200 GeV obtained dEx/dyly^o = 
620 ± 33 GeV [ 2, for the 5% most central events, our final value is lower by a 
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factor of more than two when assuming that hadronization does not change the 
local energy density. The cutoff pq — 2 GeV used in the simulations has been a 
rather conservative assumption for the initial minijets. Soft partons with px < po, 
which are produced by the nonperturbative part of nucleus-nucleus interaction and 
also contribute to the entropy production, are completely neglected. Such partons 
may stem from the color glass condensate (CGC) [s^ and should also be included 
in the transport simulation, when they are freed from the color field. In principle, 
one should combine gluons from CGC and harder gluons from minijets production 
to give a more realistic initial condition for partons after a heavy ion collision. The 
final results after parton evolution with such initial condition is thus appropriate for 
comparisons with the experimental data. We leave this interesting topic for future 
investigations. For the moment, initial conditions of gluons from CGC and minijets 
production are assumed separately for independent simulations. 

In another idea within a saturation picture of quarks and gluons in phase space, 
presented by Eskola et al. [ 4 , the authors assumed that performing the computa- 
tion of minijets at a saturation momentum, po ~ Psat, gives an estimate of the effect 
from all momentum scales, both above and below Psat- At RHIC energy they got 
Psat — 1-13 GeV for corresponding initial minijets, and estimated the final trans- 
verse energy (IEt / dy\y^Q « 660 GeV by assuming immediate thermalization and 
adiabatic expansion. Compared with the experimental data, the calculations met 
the value of the transverse energy. The smaller the cutoff po is assumed, the larger 
will be the initial giuon density. Intuitively, a denser system may achieve faster 
thermalization than a dilute system. However, immediate thermalization might be 
a naive assumption, particularly for chemical equilibration. This strongly motivates 
us to apply the developed parton cascade to inspect the timescale of thermalization 
in dependence on the cutoff parameter pq. 

In addition, Serreau and Schiff [ El have studied the role of elastic processes 
in kinetic equilibration using minijets as well as color glass condensate as initial 
conditions and found that in all cases elastic interactions alone cannot drive the 
system rapidly to kinetic equilibrium. It is interesting to see whether the situation 
changes when the inelastic gg <-> ggg processes are included. 

2. Initial Conditions 

The production of the primary partons at the very onset of a central Au Au collision 
is assumed first as a free superposition of minijets being liberated in the individual 
semihard nucleon-nucleon interactions. We introduce an additional formation time 
for every minijet, Aiy = coshy Ar^ « coshy • l/pr, which models the prior off-shell 
propagation of the to be freed, on-shell partons in the scattering process. Within 
that time span, one assumes, for simplicity, that the still virtual parton does not 
interact and moves with speed of light. The po dependence of initial conditions is 
shown in Table ^ The calculations are performed by multiplying the number of 
individual nucleon-nucleon collisions for a Au Au collision and the number of gluons 
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Table 1. Gluon number and transverse energy per unit rapidity at midrapidity 
extracted from different initial conditions and at final time 4 fm/c after cascade 
simulations. 



minijets/CGC 
Pa/Qs [GeV] 


^(init.) 


^(4Wc) 


^(init.) 
[GeV] 


^(4 fm/c) 
[GeV] 


Po = 2.0 


181 


352 


478 


272 


Pa = 1-5 


537 


663 


1082 


514 


Po = 1-4 


688 


781 


1301 


626 


Pa = 1-3 


889 


930 


1569 


775 


Qs = 1.0 


830 


528 


550 


314 



produced in a nucleon-nucleon collision at y/s = 200 GeV. The former is obtained 
by applying a Glauber picture with a Woods-Saxon nuclear distribution. The latter 
is evaluated employing the Gliick-Reya-Vogt parametrizations [E] for the parton 
structure functions. One sees that the initial gluon number and transverse energy 
increase with decreasing po . 

For the initial gluon distribution from the saturation picture we employ an 
idealized and boost-invariant form proposed by Mueller [ 7 and express it as 

f{x,p%^o = ^ - 5{p.) Q{Ql - pI) , (1) 

which is described by a bulk scale Qa^ the momentum scale at which gluon dis- 
tribution saturates. The boost-invariance leads to the equality of momentum and 
space-time rapidity, i.e., 7? = y, for the initial gluons. Distribution has been used 
in [ |S] to investigate kinetic equilibration of gluons in the subsequent evolution by 
means of a non-linear Landau equation. In contrast to the minijets the saturation 
model gives the distribution of partons with pT < Qs- We read off the value of 
parameters taken in [|H|: A^c = 3 for SU(3), c = 1.3, = 0.3, Qs = 1 GeV and the 
corresponding formation time, at which gluons become on-shell, to be r/ — 0.4 fm/c 
~ 1/Qa, and assume that gluons will be produced within a transverse radius of 6 fm 
(Au nucleus) and within \ri\ < 3 longitudinally. Then we obtain dNg/dy\y=Q — 830 
and (IEt / dy\y=o — 550 GeV, which are also included in Tabled for comparisons. 

3. Three-body Gluonic Interactions 

The three-body gluonic interactions are described by the matrix element [|5] 

'•^"^^''=(¥ (qi + m^)0 ( ki[(k/' q?)''^ + ml] ) ^^'^^^ ' ^"^^^^ ' 

(2) 

where = Anas- q± and denote, respectively, the perpendicular component of 
the momentum transfer and that of the momentum of the radiated gluon in the cm. 
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frame of the collision. We regularize the infrared divergences by using the Debye 
screening mass rnj^ which is calculated locally over the present particle density 
obtained from the simulation. The suppression of the radiation of soft gluons due 
to the Landau-Pomeranchuk-Migdal (LPM) effect, which is expressed via the step 
function in Eq. |(5J), is modeled by the consideration that the time of the emission, 
~ ^i-coshy, should be smaller than the time interval between two scatterings or 
equivalently the gluon mean free path Ag. This leads to a lower cutoff of k± and to a 
decrease of the total cross section. The dependence of the total cross section on the 
LPM-cutoff is approximately logarithmic and thus its sensitivity to the dynamics 
is only moderate. 



4. Results 

The space time propagations of gluons are simulated until a final time 4 fm/c 
at which the energy density at midrapidity at po = 2.0 GeV has already decreased 
nearly to the critical value of 1 GeV/fm'^ being characteristic for hadronization. The 
gluon number and transverse energy per rapidity at the final time are shown in Table 
n Comparing with the initial values the enhancement of gluon number is weaker 
at smaller po, which implies larger initial fugacity at smaller po- On the contrary, 
the initial gluons from the saturation model are, surprisingly, oversaturated. This 
can be interpreted as follows. According to the gluon distribution the initial 
effective temperature is T(t/) = e/3n = 2/9 * Qg- We obtain then the effective 
gluon fugacity at the formation time 

dN 

" n,, ~ 2(7V2 - l)^ - 43 asN^TfQ, ' ^""^ 

which would be valid when assuming immediate kinetic equilibrium at t = Tf. For 
the setup of parameters given before one has Ag(r/) — 8 ! So after the gluons 
are freed from the color field and become on-shell, ggg — > gg collisions will be the 
dominant processes driving the system into chemical equilibrium. 

From Tabled one also sees that in all cases of initial conditions the transverse 
energy per unit rapidity decreases with the progressing time due to the onset of 
collective expansion. The final values of dET/dy\y=o are, however, quite different. 
We see that the value at po — lA GeV would meet the experimental data. Thus the 
minijets production at pq — lA GeV seems to give an appropriate initial condition 
of gluons at RHIC. To show the decreasing of the transverse energy at midrapidity 
we depict in Fig. ^ the normalized transverse energies obtained from simulations 
with the different initial conditions. The denominators are the transverse energies 
at midrapidity at time 0.5 fm/c, which are nearly the initial values shown in Table 
n We see a unique behavior of the time evolutions of the normalized dEx I dy\y^Q^ 
especially for initial conditions employing minijets with po = 1-3 — 1.5 GeV and 
CGC with Qs — 1.0 GeV. For these initial conditions the decreasing of the trans- 
verse energy is scaled with the initial value. This indicates that the work done 
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Fig. 1. Time evolution of nomalized transverse energy per unit rapidity at midra- 
pidity. Results are obtained from simulations with different initial conditions of 
gluons, assumed in a central Au Au collision at RHIC energy. 

by the pressure is also scaled with the initial pressure when it is built up due to 
thermalization. When looking at the time evolution of collision rate per particle, 
depicted in Fig. [3 for all scattering channels, we see that the rates are in the same 
range at late times. This seems to be the reason for the scaling behavior observed 




Fig. 2. Time evolution of collision rate extracted in the central region taken as an 
expanding cylinder with a radius of 1.5 fm and within a unit interval of space time 
rapidity 77 around the collision center 77 = 0. 

for dET/dy\y=o{t), since the total collision number per unit time, which corresponds 
to the work done by the pressure, is proportional to the present local gluon density 
which is not very different from the initial value for the cases at pq = 1.3 — 1.5 
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GeV, seen in Tabled On the other hand one should realize that small angle scat- 
terings do not really contribute much to (quasi)hydrodynamical work. Transport 
cross sections, or better, transport mean free path should and will be discussed in 
detail within the parton transport approach in future investigations. 

In Fig. ^ wc also depict the normalized transverse energy per unit rapidity 
for an ideal hydrodynamical expansion by the thin solid line. Comparisons show 
that the collective expansion due to the pQCD interactions is quasi ideal at 0.5 — 2 
fm/c. At late times the decrease oi (IEt / dy\y=i^ slows down, since the collision rates 
become smaller, especially in the outer, transversally expanding region where the 
gluon density is smaller compared to the central region. 

Inspecting Fig. |21 again and focusing on the cases at po — 1-3 ~ 1-5 GeV one 
finds that, at a certain time, the smaller po is, the larger are the rates. This is 
consistent with the finding shown in Fig. 11 in [ 1 that at equilibirum the collision 
rate is proportional to temperature T which is larger in denser system initialized 
at smaller pq. The cross section is then approximately inverse proportional to . 
Therefore, at a certain time, the cross sections (^gg—^gg for gg gg and cTgg^ggg 
for gg —^ ggg scatterings are smaller at smaller pq. For instance, comparing with 
the time evolution of cross sections calculated at po — 2.0 GeV, Ogg^ggijigg^ggg) 
obtained from the simulation at po = 1-4 GeV increases from 0.5(0.3) mb at 0.5 fm/c 
to 2.5(1.4) mb at 4 fm/c respectively. The values at the final time are approximately 
1.6 times smaller than those at po = 2.0 GeV, which is consistent with the ratio of 
r2 at 4 fm/c [T = 215(262) MeV at po 2.0(1.4) GeV]. 

Figure shows kinetic and chemical equilibration observed locally in the cen- 
tral region. The kinetic equilibration is characterized by the time evolution of the 

3.5 I — ■ — , — ■ — , — ■ — , — . — , — ■ — , — ■ — , — ■ — , — ■ — I 10^^ — , — ■ — , — ■ — , — . — , — ■ — , — ■ — , — ■ — , — ^ 
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minijets p^^1 .5 GeV 
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t [fm/c] t [fm/c] 

Fig. 3. Left panel: Time evolution of the momentum anisotropy. Right panel: Time 
evolution of the effective fugacity \g = Rggg^gg / Rgg^ggg. Results are extracted 
in the central region of the expanding gluon systems simulated using the parton 
cascade with different initial conditions. 

momentum anisotropy, < > /2 < > , which maintains a value of one at equi- 
librium. For an extreme case of free streaming the final momentum distribution of 
particles in a local disc would be transversely directed, which leads to an infinite 
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value of the momentum anisotropy. Therefore initial free streaming is the reason for 
the increase of the momentum anisotropy at the beginning of the expansion, seen in 
Fig. I^lfrom the simulations with minijets initial conditions. The following bending 
over signals the onset of equilibration. We also see that full kinetic equilibrium 
comes slightly sooner, at 1 — 2 fm/c, at smaller pq. The timescalc tends to saturate 
at smaller pq. For the case using CGC as initial condition, the momenta of gluons 
are initially directed transversely in local discs according to Eq. Q). Therefore the 
momentum anisotropy has a large initial value. Surprisingly, from 0.8 fm/c the time 
evolution of the momentum anisotropy is almost identical with that obtained from 
simulations at smaller pq. 

However, we note that the maintenance of the momentum isotropy is not nec- 
essaryly equivalent to kinetic equilibrium, because even a free streaming keeps the 
momentum isotropy, but leads the system out of equilibrium. Therefore, one has 
to look at the behavior of the pT spectrum in time. While a free streaming can 
not change the pr spectrum, the spectrum steepens continuously when the system 
maintains equilibrium, since the thcrmodynamical work outwards cools down the 
system. When we compare the pr spectra at 2.0 and 4.0 fm/c, depicted in Fig. 01for 
the case when a color glass condensate is assumed as the initial condition of gluons, 
we realize that the local system at the central region is indeed at equilibrium at late 
times. One finds the same result also for minijets initial conditions as shown in more 
detail in Fig. 44 in reference [ 1 . At late times, for the transversely outer region 
the situation is different. The local system becomes out of kinetic equilibrium when 
one goes transversely away from the center to the edge of the system. That is the 
reason why the total transverse energy at midrapidity decreases only moderately at 
late times, as shown in Fig. m 

In the right panel of Fig. Owe show the ratio of the collision rates for ggg gg 
and gg —^ ggg scatterings, which is equal to the gluon fugacity Xg when the system 
is at kinetic equilibrium. The collision rates have been already shown in Fig. [3 
We realize that the timescale for full chemical equilibrium depends on the initial 
conditions. At po = 2.0 GeV the system is initially undersaturated and reaches 
full equilibrium at 3 fm/c. At smaller po the system becomes denser and chemical 
equilibrium is achieved almost initially. For a CGC situation the initial system is 
oversaturated. Full chemical equilibrium follows at a time w 1.5 fm/c. 

Finally, looking at the transverse momentum spectrum depicted in Fig. 01 we 
observe that the region around the hard scale Qs — 1.0 GeV, which is mostly 
populated by initial gluons, is less and less occupied with progressing time. The 
dominantly populated region moves to smaller pr scale as cooling proceeds. At 
2 fm/c the spectrum possesses a full thermal shape. Thermalization somehow re- 
sembles the idealistic bottom- up scenario [ 10 . However, in clear contrast to that 
scenario, the gluon number decreases for the particular parameter sets of the initial 
condition. 
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Fig. 4. Transverse momentum spectrum in the central region. 
5. Conclusion 



We have studied local thermalization of gluons in a central Au Au collision at RHIC 
energy. Initial conditions of gluons are assumed by minijets production at different 
cutoff po = 1.3 — 2.0 GeV and by a color glass condensate with the saturation scale 
Qs = 1.0 GeV respectively. We apply the newly developed parton cascade including 
99 ^ 999 to simulate the space time evolution of gluons. The numerical results 
show that kinetic equilibrium comes sooner in denser system and the timcscale 
tends to saturate, while the dependence of the chemical equilibration on the initial 
conditions is much stronger. 

Future investigations will concentrate on the energy loss of high-energy partons, 
parametrical dependence of the timescale of thermalization on Qg and concerning 
the bottom-up picture, generation of elliptic flow V2 in noncentral collisions and 
viscosity in the plasma. Moreover, quarks will be included into the parton cascade. 
A special interest will be put on the study of elliptic flow of heavy quarks. 
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